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ABSTRACT 

We study numerically the formation of molecular clouds in large-scale colliding flows including self- 
gravity. The models emphasize the competition between the effects of gravity on global and local scales 
. . . in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity - triggered by 

' a combination of strong thermal and dynamical instabilities - causes cores to form. The dynamical 

\ instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive 

i protostellar cores of a few 100 M©. The forming clouds do not reach an equilibrium state, though 

the motions within the clouds appear comparable to "virial" . The self-similar core mass distributions 
pi |. derived from models with and without self-gravity indicate that the core mass distribution is set 

O ' very early on during the cloud formation process, predominantly by a combination of thermal and 

^0 ' dynamical instabilities rather than by self-gravity. 

. Subject headings: instabilities — gravity — turbulence — methods:numerical — stars:formation — 

ISM:clouds 
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1. RAPID STAR FORMATION 

There is increasing evidence that star forma- 
tion in the solar neighborhoo d follows rapidly upon 
molecular cloud formation (|Hartmann et al.l 1200 It 
iBallesteros-Paredes fc Hartmamj |2007| and references 
therein). This evidence suggests that the density en- 
hancements in which stars form are produced during the 
cloud formation phase; thus understanding cloud forma- 
tion is essential to understanding star formation. More- 
over, it appears that non-linear density perturbations 
need to arise quite early in cloud formation, as mas- 
sive, finite molecular clouds are highly susceptible to 
large-scale gravitational collapse which could overwhelm 
small, stellar-mass fragmentation (jBurkert fc HartmannI 
|2004 ). While several investigations have adopted var- 
ious assumed forms of initial and/or driven turbulent 
motions to produce the n ecessary small-scal e structure 
(iKlessen fc BurkertI 120001 : jPado a n fc Nordlmldl l2002t 
iBate et al.l l2002l I2003D . it is preferable to have these 
structures arise naturally. Thus, a close look at insta- 
bilities in cloud formation which could lead to strong 
densi ty fluctuations is neede d. 

iBallesteros-Paredes et a l.l ()1999aD and lHartmann et ahl 
()2001h proposed that cloud formation as the result of 
pileup of material by large-scale flows is an essential 
mechanism for explaining the "crossing time problem" , 
i.e. the observation that the typical age spreads in the 
stellar populations of many large star-forming regions 
are often substantially smaller than the lateral crossing 
timescales; in the large-scale flow picture, no information 
is transmitted laterally, i.e. perpendicular to the large- 
scale flow. This picture works only if star formation fol- 

^ Dept. of Astronomy, University of Michigan, 500 Church St., 
Ann Arbor, MI 48109-1042, U.S.A 

^ Oxford University, Astrophysics, Denys Wilkinson Building, 
Keble Road, Oxford, OXl 3RH, United Kingdom 

^ Universite Claude Bernard Lyon 1, CRAL, Observatoire de 
Lyon, 9 Avenue Charles Andre, 69561 St-Genis Laval Cedex, 
France; CNRS, UMR 5574; ENS Lyon 

Universitats-Sternwarte Miinchen, Scheinerstr. 1, 81679 
Miinchen, Germany 



lows closely upon molecular gas formation. Expanding H 
II regions, supernova bubbles, and spiral density waves 
are all obvious candidates for large-scale supersonic flows 
which can sweep up material, and there is considerable 
direct observational evidence for rapid sta r formation in 
these environments (jHartmann et al.ll2001l and references 
therein). Thus a plausible place to look for stellar core- 
forming instabilities is in the post-shock material of the 
large-scale flows. 

Several numerical studies relevant to this problem have 
now be en undertaken (see the discussions of the litera- 
ture in flleitsch et al.ll200"6l and 'Vazq uez-Semadeni et al.l 
I2007D . These calculations typically assume converging 
flows to keep the shocked gas within the computational 
volume, but this can easily be extended to describe a 
more generic situation by recasting the problem in the 
rest frame of the shock(s). These models overcome the 
limitations of p revious turbulent fragmen tation models 
(see review by iMac Low fc KlessenI |2004[ ) by avoiding 
ad hoc assumptions about the source of turbulence and 
boundary conditions. Indeed, the converging flow mod- 
els demonstrate with ease that flows provide a natural 
mechanism for the generation of structure and turbu- 
lence in clouds (Audit & HcnnebcUe 2005: Heitsch et alj 
2005: Vazauez-Semadc ni ct al. 2006; Heitsch ctaL 200g 
Hennebelle et al.ll2007l : iHennebelle'fc'AudidlM)7h . 

Although this study is motivated by the sce- 
nario of molecular clouds be ing transient entities 
(jBallesteros-Paredes et allll999a| ). forming and dispers- 
ing in background flows within a few free-fall times, 
the results presented here are not restricted to this sce- 
nario. Colliding flows can appear even in large scale 
gravitational instabilities, linking our models to the 
alternative scenario of "Giant Molecular Clouds" liv- 
ing for substantially longer than only a few free-fall 
times. The issue of cloud lifetimes is currently a mat- 
ter o f debate (e.g. iKrumholz fc TanI 120071 : lElmegreenI 
120071 iMcKee fc Ostrikeij 120071 ), and s eems to depend 
strongly o n the galactic environment (jHartmann et al.l 
2001; Heit sch et al.ll2006h . For reasons discussed in |5l 
our models cannot predict cloud life times. Thus, the 
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emphasis of this study is on the onset of star formation. 

The purpose of this study is to compare the fragmen- 
tation processes in simulated converging flows with and 
without self-gravity in order to show that the rapid on- 
set of star formation is pretty much unavoidable within 
the scenario where molecular clouds form i n con verg- 
ing flows. Recently. IVazquez-Semadeni et al.l ()2007D pre- 
sented a study of star formation in clouds formed by col- 
liding flows, emphasizing the long-term evolution of the 
system. Here we focus more on the initial development of 
the cloud, discussing the consequences of cooling, cloud 
geometry, and gravity for the star formation process. 

Cooling and gravity both are fragmentation agents, 
with the difference that gravity can be relevant on all 
scales which surpass the Jeans length, while th e iso- 
baric condensation mode of the thermal instability (|Fieldl 
|1965() is limited to scales set by the sound speed and 
the cooling time, Ac = Cs Tc (e.g. iBurkert fc Lii]|[2000l : 
iHennebelle fc Audit! I2007f) . Thus, in the early stages of 
cloud formation, when only a little mass has accumu- 
lated, the thermal instability is the dominant fragmen- 
tation agent. 

Non-linear density pertu rbations collapse at a highe r 
rate than the global cloud (jBurkert fc Hartmannll200l ). 
thus allowing stars to form locally before the whole-sale 
collapse of the cloud. We find that the strong dynami- 
cal and thermal instabilities generate non-linear density 
perturbations for rapid local collapse, while still allowing 
for the global build-up of the cloud. 

The physics and methods are summarized in Sj2l fol- 
lowed by the model results in !J31 Readers solely inter- 
ested in the discussion of the results and their conse- 
quences should directly proceed to Sjll^nd SjSl 

2. PHYSICS AND METHODS 

Our study focuses on the effects of global versus local 
gravity on the one hand, and on the rapid generation of 
substructure in the colliding flows on the other. Since 
we are interested in global gravitational effects, we can- 
not use the periodic boundary cond itions of earlier tur- 
bulent fragmentatio n studies (e.g. iKlessen et al.l [20001: 
Heitsch ct al. 200 1; iPadoan et alj 120011: iGammie et all 
2003t IVazquez-Semadeni et al.ll2005D . Neither will we 
generate or drive turbulence by imposing a randomly 
chosen velocity or density field, but instead, we will rely 
on turbulence generated by the dynamical instabilities 
arising from the collisions of the fiows. 

2.1. The Models 

We ran four models, whose parameters are listed in 
Table [TJ All models are run on a fixed grid with the in- 
streaming gas fiowing along the cc-direction, entering the 
domain at the {y, z)-planes. To trigger the fragmentation 
of the (otherwise plane-parallel) interaction region, we 
perturb the collision interface. We chose the perturba- 
tions of the collision interface from a random distribution 
of amplitudes in Fourier space with a top hat distribution 
restricted between wave numbers k = 1..4. 

Model Gs (for "gravity in shell") can be interpreted 
as colliding cont inuous gas streams in spiral shocks (e.g. 
Tilanus fc Allen 1990 for observational evidence, and 
Dobbs fc BonncU 2007 for numerical models), or as a 
close-up view of two expanding and colliding super-shells 



Table 1. Model Parameters 



Name 




LxLyLz [pc] 


gravity t^„a [Myr] 


V [pc] 


Hfl 


256 X 512^ 


22 X 44^ 


no 14.5 


2.2 


Gfl 


256 X 512^ 


22 X 44^ 


yes 14.5 


2.2 


Gf2 


256 X 512^ 


22 X 442 


yes 14.5 


4.4 


Gs 


2563 


443 


yes 16.0 


2.2 


Note. 


— 1st column: 


Model name. 


2nd column: resolution 


3rd 



column: physical grid size. 4th column: gravity. 5th column: end 
time of run. 6th column: amplitude of interface displacement. 



in the LMC. The collision interface is plane-parallel, ex- 
cept for the imposed perturbations. Material is free to 
leave the box in the lateral (i.e. perpendicular to the 
inflow) directions. 

In models Gfl and Gf2 (for "gravity in finite cloud"), 
we restrict the inflow to a cylinder of elliptical cross sec- 
tion with an ellipticity of 3.3 and a major axis of 80% 
of the (transverse) box size, mimicking two colliding gas 
streams in a more general geometry. Again, the collision 
interface is perturbed. The motivation here is to gener- 
ate one finite cloud in order to study global gravitational 
effects. 

Finally, model Hfl is a non-gravitating version of Gfl, 
to compare the role of gravity versus that of the thermal 
instability for the fragmentation of the gas streams. 

The inflow density in all models is no = 3 cm~^ at a 
temperature of To = 1800 K and an inflow velocity of 
7.9 km s"^, corresponding to a Mach number oi A4 = 
1.5. The flows are initially in thermal equilibrium. The 
models start at time t = with the collision of the two 
flows. For models with spatially constrained inflows (Hfl, 
Gfl, Gf2), the fluid is at rest everywhere except in the 
colliding cylinders. 

The finite cloud models Gfl, Gf2 and Hfl have a grid 
cell size of AL — 8.6 x 10~^ pc, while model Gs has 
one of AL = 1.7 x 10"^ pc (Tab. [l]). We note that this 
does not constitute the physical resolution power of the 
simulation. At minimum, the stencils (i.e. the support 
points) used for the higher-order reconstruction of the 
fluid states at the cell wall will render the cells within 
one stencil not independent. In other words, conclusions 
should not be drawn from structures of 4 cells or less of 
linear size. Thus, we use only cores with 64 or more cells 
for analysis. 

2.2. Boundary Conditions 

The x-boundaries are partly or entirely defined as 
inflow-boundaries, depending on the model. Indeed, the 
inflow is either defined over the whole (y, z)-plane (model 
Gs), or within an elliptical surface (models Hfl, Gfl, 
Gf2: see 2.11 for details). The y and z boundaries, - 
as well as the part of the x-boundaries that is not oc- 
cupied by the inflow in models Hfl, Gfl and Gf2 -, are 
open, meaning material is free to leave the simulation do- 
main through these boundaries. This is bound to cause 
trouble once material tries to "come back" during the 
later stages of the simulation (note that this material is 
not actually coming back, but that it is the result of the 
extrapolation of the last active cells properties, i.e. mate- 
rial with the properties of the last active cell layer within 
the domain will try to enter the domain) . This inevitably 
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will happen once global gravity dominates over the over- 
pressurized material shooting out of the flow-collision re- 
gion. However, the "re-entering" material does not reach 
the central cloud region within the simulation time, and 
in any case contributes only a negligible amount to the 
total mass within the box. 

The situation becomes more critical once material is 
leaving the simulation domain in the x-direction, i.e. 
once it is moving against the inflow. Since the bound- 
ing shocks (and the cooling) will set the density and the 
temperature of the post-shock gas, the physical state of 
the gas will be undefined once the bounding shocks move 
off the grid. This will render the "returning" material es- 
sentially in a hydrodynamically inconsistent state. Thus, 
once material encounters the cc-boundaries, the simula- 
tion needs to be stopped. 

In model Gs, the inflow velocities are slightly reduced 
at the edges of the domain, mimicking the velocity proflle 
of an expanding shell of material driven by two sources at 
a distance of approximately 100 pc to the left and to the 
right of the mid plane of the simulation domain. This 
reduces the amount of material collected at the edges 
of the domain, and thus limits the edge effects due to 
gravity at later stages of the simulation. 

2.3. Hydrodynamics and Atomic Line Coolants 

As in our previous studies of colliding flows, we 
used the higher-order gas-kinetic grid method Pro- 
teus (jP rcndcrgast & Xu 1993; Slvz & Prcndcrgast 
Heitsch et al.ll2004t iSlvz et al.l l2005: Heit sch et al] 



1999; 



2006 



20071) . allowing full control of viscosity and heat con 
duction. The code evolves the Navier-Stokes equations 
in their conservative form to second order in time and 
space. The hydrodynamical quantities are updated in 
time unsplit form. 

The heating and cooling rates are restricted to op - 
tically thin atomic lines following IWolfire et al.l (|1995f ). 
Dust extinction becomes important above column den- 
sities of A^(HI) « 1.2 X lO^^cm"^, which are only 
reached in the densest regions modeled. Thus, we use 
the unattenuated U V radiation field for grain heating 
(jWolfire et al.lll995l ). expecting substantial uncertainties 
in cooling rates only for the densest regions. The ioniza- 
tion degree is derived from a balance between ionization 
by cosmic rays and recombination, assuming that Ly a 
photons are directly reabsorbed. Numerically, heating 
and cooling is implemented iteratively as a source term 
for the internal energy e of the form 



dfc = nT{T) ~ n^K{T) [erg cm'^ s' 



(1) 



Here, F is the heating contribution (mainly photo- 
electric heating from grains), nA the cooling contribution 
(mainly due to the CH HFS line at 158/Ltm). Since the 
cooling and heating prescription has to be added outside 
the fiux computations, it lowers the time order of the 
scheme. To speed up the calculations, equation ([T]) is 
tabulated on a 2048^ grid in density and temperature. 
For each cell and iteration, the actual energy change is 
then bi-linearly interpolated from this grid. 

The cores forming due to gravitational collapse reach 
densities of a few 10^ cm~^, far beyond the applicable 
range of our cooling curve. Strictly speaking, we should 
therefore extend the cooling curve to include molecu- 
lar lines at high densities. However, our cooling curve 



reaches an equilibrium temperature of approx 12K for 
n > 10'^ cm~'^, close enough to a realistic temperature 
for molecular cores. Since we cannot resolve the core 
structure anyway, we chose to stick to this simplified 
treatment. 

The cooling curve is limited to densities of n < 
10^ cm~'^, to prevent a "catastrophic" collapse which 
would be generated by a sub-isothermal effective equa- 
tion of state. The sudden reversal to an adiabatic equa- 
tion of state helps limit the densities and prevents nu- 
merical artifacts caused by single, very high density cells. 
However, this stiffening of the equation of state - if intro- 
duced at too low densities - could stabilize the cores, pre- 
vent their fragmentation and render them more prone to 
dispersion. We experimented with the density threshold 
n-max and found that a value of Umax — 10^ cm^'^ pre- 
vents the run-away collapse while allowing the cores to 
remain small (and dense) enough to stay gravitationally 
bound once they have formed. See §3 3.41 for a discussion 
of the resolution limits. 

2.4. Gravity 

Self-gravity is implemented as an external source term 
in time-unsplit form. The Poisson-equation is solved 
via a non-periodic Fourier solver, using the (MPI- 
parallelized) fftw (Fastest Fourier Transform in the 
West) libraries. We tested this against direct summa- 
tion to assure that the Poisson equation is solved accu- 
rately. The non-periodic solver needs twice the grid size 
for padding the Fourier transforms. This limits the res- 
olution of our simulations to effectively 512'^ cells. 

2.5. Core Identification 

We use two methods to identify cores in the model 
data. To find gravitationally b ound objects, we emp loy 
the CLUMPFIND algorithm (William's e^aD [THj in 
a modified version (Klcssen et al. 2000). We then test 
whether the structures identified by CLUMPFIND are 
gravitationally bound or collapsing, by checking their 
Jeans mass, the ratio of thermal and (internal) kinetic 
energy over gravitational energy, and the velocity diver- 
gence. If all three tests are passed, a structure is accepted 
as a core. We track individual cores by identifying the 
closest "neighbor" to a given core in the next timestep 
(where "timestep" does not mean the CFL-timestep, but 
the time between writing data sets). The simulations 
presented here form sufficiently few cores for this simple 
method to be accurate. 

The second method is a simple clipping algorithm, 
motivated by the fact that due to the thermal insta- 
bility, dense coherent regions are generally well defined 
in our models. The method selects the maximum den- 
sity and builds a tree structure around the central cell, 
thus connecting all cells above the given density thresh- 
old of utii > 50 cm^"^. Once this threshold is reached, 
the process restarts with the next-lower density peak not 
included in the previous structure. The resulting cores 
are accepted independently of whether they are gravita- 
tionally bound or collapsing. 

3. MODEL RESULTS 

The general signature of fast local fragmentation in col- 
liding flows is most easily recognized in the morphologies 
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of the clouds ( ^3 3.1[) . A more quantitative measure can 
be gleaned from the core mass evolution and the energy 
distribution ( H'S 3.2p . Dynamical signatures are discussed 
in ^3 3.31 and i )3 3.41 sounds two cautionary notes regard- 
ing the numerical resolution. 

3.1. Morphologies 

We begin by comparing the morphologies of the clouds 
forming in colliding flows (models Hfl, Gfl and Gf2, 
§3 3.1 3.1.ip . A new formation mechanism for massive 
cores is discussed in H3 3.1 3.1.21 . 

3.1.1. Global Collapse and Filament Formation 

The top row of Figure [1] shows three time instances 
of model Hfl, seen along the inflow direction. The col- 
liding flows cause a big "splash" , the effects of which 
are still noticeable 7.6 Myr after the initial flow colli- 
sion (left column), however, rapid cooling leads to strong 
density enhancements in the interaction zone, and in 
combination with the dynamical instabilities triggered 
by the perturbed interface, to strong fragmentation. 
Note the radial fllaments and the outermost "bound- 
ing ring" at 7.6 Myr. The radial fllaments are similar 
to tho se seen in the co nceptually equivalent models by 
IVazgu ez-Scm adeni et al.. (2007): the compressed mate- 
rial is escaping the interaction region by the way of least 
resistance, i.e. laterally to the inflow. Once an "escape 
channel" has been formed due to small fluctuations in 
the external pressure, the resulting pressure deflcit will 
ensure that the channel will continue to be used by subse- 
quent material. The ring-like structure is just the shock 
wave from the initial splash caused by the colliding flows. 

With increasing time, more and more mass is collected 
in the interaction region, rapidly assembling the cloud. 
Note that for the sake of greater detail the center and 
right-hand column of Figure [T] show only 3/4 (in linear 
extent) of the simulation domain corresponding to a box 
length of 33 pc. The non-gravitating model Hfl also 
continues to collect mass as time proceeds, however the 
column densities reached do not exceed a few 10^^ cm~^. 

This changes with the introduction of gravity (center 
and lower row of Fig. [T]). At 7.6 Myr, the structures are 
still pretty much indiscernible, while at 11.4 Myr, the 
flrst regions of high column density have formed. While 
the cores forming in the cloud result from local collapse 
(see below), the filaments at later stages (14.5 Myr) are 
a consequence of the global collapse of the whole cloud. 
This can be more easily seen in the right-hand column of 
Figure m The same models at the same times as in Fig- 
ure [T] are shown, but now seen perpendicularly to the 
inflow. Clearly, at late times, the initially elongated 
"red" structure (the dense part of the cloud) crumples 
under its own weight. Because of the (generalized) non- 
circular inflow cross section, this leads to the formation 
of a fllament. This mechanism for filament formation 
offers a substantial reservoir of mass for further star for- 
mation. The subsequent fragmentation of the filaments 
is possibly enhanced by t he thermal instability (see e.g. 
iTsuribe fc Inutsukal[200l . 

The larger amplitude of the collision interface pertur- 
bation in model Gf2 (see Tab[T]) leads to a stronger initial 
fragmentation, mirrored in a more distributed core for- 
mation at later stages. 



Note that the high-density regions (in Fig. [T]) do not 
necessarily form at the center of the cloud. On the 
contrary, there is a tendency for material to collect 
away from the center, forming filaments (model Gf2 
at 14.5 Myr), or at least extended dense cores. This 
is a mild version of t he edge effect in collapsing finite 
sheets, as discussed bv lBurkert fc HartmannI (|2004l ) and 
applied to a model of the O rion star forming region by 
iHartmann fc BurkertI (|2007D . Figure [3] is more specific 
about the actual mechanism: It shows a map of the (pro- 
jected) gravitational accelerations |V$| in the midplane 
perpendicular to the inflow (to be compared to the bot- 
tom row of Fig. [T|) . Contours denote the column density, 
and the actual potential gradient — V$ is indicated by 
the arrows. The color table is identical to Figure [l] i.e. 
yellow/red indicates strong accelerations. The ring-like 
structure of strong accelerations towards the center is ob- 
vious. Note that the accelerations extend towards larger 
radii than their accompanying high-density structures: 
material is accelerated at the edges and is being piled up 
further down the (radi al) flow. 

The simulations of IVazguez-Semadeni et al.l ()2007| ) 
show a similar effect as the models discussed here, 
although much stronger. This quantitative differ- 
ence might be a consequence of the choice of ini- 
tial perturbations in th e two simulation sets. While 
I Vazquez- Semadeni et al.l (|2007l ) put small-scale pertur- 
bations in their inflow velocities, we perturb the collision 
interface, which in turn leads to a stronger excitation 
of dynamical instabilities. The predominant instability 
arising in our setup is t he non-linear thin shell instability 
(NTSI. IVishniadll994f ). in combination with shear flow 
instabilities and the thermal instability (|Fieldlll965l ) , pre- 
venting the formation of a more or less uniform slab. 

3.1.2. Dynamical Focusing and Massive Core 
Formation 

The NTSI is driven by the ram pressure imbalance at 
the troughs of the rippled interaction surface: the con- 
cave side will have an excess of ram pressure because the 
flows are focused into the troughs, while the convex side 
will experience a deficit of ram pressure because the in- 
coming gas is deflected. Thus, the NTSI pr ovides a very 
efficient mechanism to collect gas locally. iHueckstaedtl 
(|2003l ) demonstrated this effect with the help of self- 
gravitating two-dimensional models. 

This focusing effect still holds in three dimensions, as 
can be seen from the time sequence of logarithmic col- 
umn density maps of two colliding flows under the effect 
of self-gravity, as shown in Figure ID The leftmost panel 
at t = 0.8 Myr shows structures still pretty close to the 
initial conditions. The cooling has not yet led to percep- 
tible fragmentation, and turbulence has not yet devel- 
oped. At 8.0 Myr, the NTSI is already in full swing, and 
on the left edge of the cloud, roughly in the mid-plane, 
the first core starts to form. This core is located at one of 
the troughs amplified by the NTSI, so that it had ample 
opportunity to collect instreaming material. This seems 
to continue all the way up to 16 Myr, at which point the 
cloud has grown globally unstable (see 3.2p . indicated 
by the frenzy of core formation all over the cloud. We 
will present a more detailed discussion of the dynami- 
cal focusing in a subsequent paper, including models at 
higher resolution. 





Fig. 1. — Time sequence of logarithmic column density maps for models Hfl (top), Gfl (center) and Gf2 (bottom), seen along the inflow 
direction. At t = 7.6 Myr (left column) the full domain (44 pc) is shown, while we restrict the field of view to the central 3/4 (33 pc) of 
the domain at later times, to highlight the small-scale dense structures forming. 



Isolating the gravitationally collapsing cores with 
CLUMPFIND (see ^223| and tracking the core masses 
with time yields Figure El The dynamical focusing seems 
to be a very efficient mechanism to collect substantial 
mass in a small volume over a short time: 9 Myr after 
initial flow contact, the first 100 MQ-core has formed. At 
later times, the mass accretion rates get steeper: there 
is more mass available, and (concurrently) the potential 
well deepens, so that more massive cores can form within 
shorter times. The first core does not participate in the 
steepened accretion history because it sits at the edge of 
the potential well so that it does not benefit from the 



higher densities. 

Obviously, stellar feedback will influence the cloud dy- 
namics and star formation efficiency at the late stages of 
the cloud evolution. 

3.2. Collapse History 

3.2.1. Energy Equipartition 

A global measure for the cloud evolution under the 
effect of gravity is the equipartition parameter 
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Fig. 2. — Time sequence of logarithmic column density maps 
for models Hfl (top), Gfl (center) and Gf2 (bottom), seen perpen- 
dicular to the inflow direction. The full computational domain is 
shown, measuring 22 X 44 pc. 



i.e. the ratio of the total kinetic and thermal energy over 
half the total potential energy^. We neglect any surface 
terms in the determination of aeq, which only allows us 
to obtain a rough approximation of the "true" energetic 
state of the cloud. Nevertheless, as a closer study of Fig- 

^ We decided to follow IBallesteros-ParedesI I I2006I ) and replace 
the term "virial parameter" , since this usually has the connotation 
of "virial equilibrium" , an assumption, which generally holds only 
for an ensemble of molecular clouds , or for a tim e-average of one 
cloud over many dynamical times ijMcKeel 119991 ) . Since we are 
concerned here with a single molecular cloud on short timescales, 
the notion of virial equilibriu m is inapplicable — while th at of energy 
equipartition may still hold l lBallesteros-Paredesll2006l) . 



ure [6] demonstrates, the time evolution of aeq mirrors 
the cloud morphologies (Figs. [4] and [T]). Figure [6] shows 
for all the gas above a given threshold density. A 
low threshold density means that most of the mass and 
most of the volume enter the calculation of Oeq. Increas- 
ing the density threshold emphasizes more and more the 
dense cores that form later. The lowest density thresh- 
old is set to nth — 10^ cm~^, since we are interested in 
the equipartition parameter of the isolated cloud, and 
not of the total simulation volume including the (highly 
energetic) inflows. The following items are noteworthy: 

(1) The dense regions tend to be gravitationally bound 
{aeq < 1), while the global cloud behavior (illustrated 
by the lowest-density curves) approaches aeq ~ 1 by the 
end of the simulations. Thus, in terms of observables, the 
cloud exhibits an aeq consistent with "virial equilibrium" 
to within a factor of two, well within observational uncer- 
tainties in terms of mass estimates, velocity dispersions 
derived only from line- of-sight motions, and the ehmina- 
tion o f surface terms f Ballesteros-Paredes et all Il999bl : 
Balle steros-Pare des 200(f^ 

(2) The solid lines indicating the lowest density thresh- 
old of nth = 10^ cm~^ reach aeq = 1 between 10 and 
13 Myr. A lower density threshold can be interpreted as 
tracing a larger volume, thus the evolution of aeq for low 
density thresholds indicates that global collapse lags be- 
hind the local collapse - isolated dense cores form before 
the cloud can collapse globally. However, the cloud does 
show the onset of global collapse. 

(3) At high density thresholds, aeq < 1 for all times (at 
which high densities are available) : the massive cores (see 
the corresponding mass history) are fully gravitationally 
unstable. Note that aeq does not drop further, although 
the mass increases: the cut-off of the cooling curve leads 
to a stabilization of the cores and prevents catastrophic 
collapse. 

(4) The clouds do not go through an "equilibrium" 
stage, but start to collapse locally during their formation: 
stars can form locally without a global collapse of the 
cloud. 

Figure [6] allows only an indirect statement about the 
scale- wise evolution of aeq- A more accurate measure is 
the ratio of the respective Fourier spectra of the kinetic 
and internal energy over the gravitational energy. This 
yields aeq{L), a scale-dependent measure of the cloud's 
stability against gravity (Fig. [7]). As in Figure [6l we 
show the three gravitational models, but now at various 
times. Note that we do not select for gas in the cloud 
or gas above a density threshold. This is because mask- 
ing leads to additional structure, which in turn results 
in "noise" signals in the Fourier spectra. To facilitate 
an easier comparison to spatial scales we plot the scale- 
dependent aeq against physical length scale, rather than 
against wave number. Since aeq is only a rough mea- 
sure of the system's energetics, the following discussion 
should be seen as a qualitative analysis. 

At early times, the system is gravitationally stable on 
all scales. This is not surprising since we are now looking 
pretty much at the whole box (we removed the L = 44 pc 
mode, since this would just be the mean). The small 
scales definitely collapse first: they are the first to fall 
beneath aeq = 1 with increasing time. The minimum 
in aeq around L = 1 . . . 3 pc at earlier times stems from 
a "conspiracy" between the kinetic and potential energy 
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Fig. 3. — Time sequence of accelerations for model Gf2, projected in plane perpendicular to the inflow. Colors denote | V'I'I, the contours 
represent density, and the arrows indicate V4>. This should be compared to the bottom row of Figure [T] 





Fig. 5. — Core mass evolution for model Gs. The global collapse 
of the cloud is paralleled by a frenzy of core formation at late times. 




Fig. 4. — Time sequence of logarithmic column density maps for 
model Gs, seen perpendicular to the inflow direction. The domain 
measures 44 pc in linear extent. 



scales: On larger scales {L >3 pc), the kinetic energy of 
the large-scale inflow still dominates the energy budget 
at early times, whereas on the small scales, the potential 
energy drops faster with decreasing scale than the kinetic 
energy. The small structures forming due to thermal and 
dynamical fragmentation did not have enough time yet 
to collect a significant amount of mass. 

In summary, we note that in all models local collapse 
wins over global collapse, i.e. the small scales gener- 
ated by thermal and dynamical fragmentation collapse 
first. Global collapse however occurs and feeds more ma- 
terial into the already active "star forming" region. The 
clouds do not reach an "equilibrium stage" , but proceed 
from formation directly to local collapse and only then 
to global collapse. This result does not depend on the 



presence or absence of stellar feedback, in contrast to the 
late-stage evolution of our model clouds, where feedback 
will have a deciding influence. 

3.2.2. Total and Core Mass Evolution 

We have already seen the core mass history of model 
Gs (Fig. [5]). Figure [5] shows the total mass evolution 
of all models. Thick lines refer to gas at T < 100 K, 
which can be identified as cloud gas due to the thermal 
instability. Thin lines denote the gas at T > 100 K. The 
symbols stand for the total mass within collapsing cores 
for each model, i.e. for Af,, the mass that constitutes 
the reservoir for star formation. 

Note that the mass is scaled logarithmically, while the 
(cold) cloud mass evolves linearly with time. Essentially, 
despite all the su bstructure, the collid ing flows form one 
cold slab (see also lHeitsch et al.ll2006l ). Comparing mod- 
els Hfl, Gfl and Gf2, gravity does not affect the global 
cloud mass: the inflows determine the mass in the cold 
gas phase (i.e. in the cloud), and thus the material avail- 
able for making stars. Model Gs differs from all others 
just because of the larger extent of the collision site. 0th- 
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Fig. 6. — The equipartition parameter Oeq against time for the three gravity models Gs, Gfl and Gf2 as indicated in the panels. The top 
half of each panel gives agq{n > n^j,), the bottom half shows the total mass M(n > nj^). The line styles stand for the threshold densities 
rifii as indicated in the top panels. 
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Fig. 7. — The equipartition parameter Oeq against scale, for the three gravity models Gs, Gfl and Gf2 as indicated in the panels. The 
line styles stand for the times at which a^q has been measured. All models have been cut at a scale of 0.3 pc, although the numerical 
resolution is a factor of approximately 3 (6) higher for model Gs (Gfl, Gf2). 
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Fig. 8. — Mass history of all models. Thick lines denote 
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styles represent the four models. Symbols refer to t he to tal mass 
in collapsing cores identified by CLUMPFIND (see ^2 2.5|l . 



erwise, its evolution regarding the cloud mass is qualita- 
tively similar. 



The core masses evolve in parts nearly exponentially. 
This is a consequence not of the individual mass accretion 
events, as Figure [5] demonstrates: those are fairly linear 
with time. Rather, the non-linear evolution is a conse- 
quence of the explosion of star formation activity once 
sufficient mass has accumulated. Comparing the mass in 
the cores, M*, against the cloud mass Md - identified 
with M{T < lOOK) we note that over a large stretch 
of time, the "star formation efhciency" M^,/Mci < 0.1 - 
the collapse is initially truly local. 

3.2.3. Core Mass Distribution 

Since the thermal instabilities are the first fragmen- 
tation agent in the colliding gas streams, will they set 
the core mass distribution in molecular clouds? Fig- 
ure [9] shows the core mass distribution in logarithmic 
mass intervals for all four models. The (fit) slopes 
dlnN/d{\nM) = s are indicated in the legend. The 
Salpeter IMF would have an exponent of s = —1.35. 
While the slopes are approximat ely consistent with 
observed core mass distriljutions ([Kramer et al.l Il996t 
ISchneider et"al][200l . we note that there is at most one 
decade in mass for the fits. Moreover, these mass spec- 
tra are most likely affected by the stiffening of the equa- 
tion of state for n > 10^ cm~'^ ( H2 2.3p . i.e. low-mass 
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objects may be under-represented. This may cause a 
flattening of the spectra of the gravity-models (Gfl, Gf2 
and Gs) - possibly visible in the lower right panel (13 
Myr) of Figure [9] -, whereas model Hfl shows a rather 
well-defined power law down to the lowest masses at late 
times. Note that the cores contributing to Figure [9] are 
not necessarily self-gravitating. Over a period of 4 Myr, 
the spectra roughly keep their shape, another indication 
that the fragmentation is mainly due to thermal effects 
rather than gravity. 

3.3. Gas Dynamics 

Figure [To] summarizes the various line-of-sight velocity 
dispersions (LOSVD) for the cold (T < 100 K) gas, for aU 
models. All LOSVDs are one-dimensional and density- 
weighted, e.g. the total LOSVD is 



3 / ndV 



1/2 



(3) 



The top panel shows the velocity dispersion along the 
inflow direction {(Jx), the transverse velocity dispersion 
Gyz, and the total dispersion. For all models, Ux is high- 
est, initially increasing slightly with time, and later lev- 
eling off. The initial rise is due to the still acting NTSI, 
i.e. the ripples in the interaction interface are still am- 
plified, while at later times, the NTSI is saturated, and 
global gravity takes over. The latter can be seen when 
comparing model Hfl and Gfl , which are identical except 
that Gfl has self-gravity. For both, ax initially evolves 
similarly, up to t w 7.5 Myr. After that, fTa;(Gfl) starts 
to level off, while ax(iiil) still continues to rise: gravity 
constrains the slab and suppresses the further growth of 
the NTSI. This is obvious in Figure El 

The transverse velocity dispersion ayz (medium lines 
in top panel of Fig. [TO)) shows the influence of global 
gravity at later times for models Gfl, Gf2 and Gs, while 
that of model Hfl drops with time. For model Hfl, the 
total energy input of the inflows is balanced by radiative 
losses and out-streaming material: the total velocity dis- 
persion stays approximately constant with time. Not so 
for the self-gravitating models: they all have increasing 
total velocity dispersions with time. 

In the bottom panel of Figure [TUl the thick lines stand 
for the total velocity dispersion again (as in the top 
panel). The thin lines {{(Jcore)) denote the velocity dis- 
persion within each core, averaged over all cores. This is 
the same quant ity as the "irL t ernal" velocity dispersion 
as discussed by iHeitsch et al.l ()2006[) . Due to the strong 
radiative losses, this velocity dispersion is subsonic (the 
sound speed in the cold gas is Cg ~ 0.3 km s^^). In 
contrast, the velocity dispersion of all gas in the cores, 
(Tcore IS slightly supcrsouic, and increases with time due 
to global gravitational collapse. 

3.4. Two Comments on Resolution 

The left column of Figure [Tl] shows the histograms 
of the Jeans length for all three self-gr avitating mod- 
els. T he dashed vertical line indicates the lTruelove et all 
()1997() criterion with a safety factor of 4, mandating a 
minimum number of cells per Jeans length. We evaluate 
the criterion locally, i.e. per cell. Clearly, with increas- 
ing time, more and more cells fall below the resolution 
limit (to the left of the vertical dashed line). A more 



detailed view is offered by the set of panels on the right 
side of Figure [TlJ They show the Jeans length in each 
cell against the corresponding density. Again, dashed 
lines denote the Truelove limit, the upper line for 4 cells, 
the lower one for 2. Although only a minor fraction of 
cells is unresolved, they do exist. The strict correlation 
between n and Aj for n < 10^ cm~^, and the scatter 
at larger densities, is a direct consequence of the cooling 
curve: we are far up the isothermal branch, so that the 
thermal timescales are much shorter than the dynamical 
timescales. As discussed above, we switch off the cooling 
for densities n > 10^ cm~'^, to prevent unphysically high 
densities and subsequent numerical problems. The effec- 
tive equation of state reverts to adiabatic with 7 = 5/3 
at that point, so that the temperature increases due to 
the strong compressions. 

There are several reso lution criteria for thermall y un- 
stable systems (see e.g. iHennebelle fc AuditI l2007l for a 
discussion). The dynamically most stringent condition 
is to resolve the cooling length Ac = CsTc. If this length 
scale is not resolved in the cold gas, the (isobaric) con- 
densation mo de of the the rm al instability wifl be u nder- 
estimated (see lFieldlfl965l and iBurkert fc Linll2000t l . The 
bulk of the cold gas resi des at 40 K for our cooling curve 
(see 'Heits ch et aLll2006[ ). and the cooling time scales in 
the cold gas are on the order of 10"* years, so that the 
cooling length scale is Ac ~ 6 x 10~^ pc, slightly be- 
neath the nominal resolution for models Hfl, Gfl and Gf2 
(see §2 2.ip . Underestimating the condensation mode of 
the thermal instability will result in fewer low-mass frag- 
ments. In that sense, the core mass budgets ( W3 3.2p and 
the fragmentation history in our models are conserva- 
tive estimates: at higher resolution, more fragmentation 

- and possibly earlier low-mass star formation - is ex- 
pected. 

4. DISCUSSION 

4.1. Global Collapse and Filament Formation 

Galactic star formation seems to prefer filamentary 

- at le ast elongated - ra ther than sph erical molecular 
clouds (ISchneider fc Elm cgrccn 1979; Ch ini et all 119971: 
Hartmann et al.ll2001l: " lHatcheU et al. 2005nAlves et all 
200711 . The finite extent of the forming cloud opens two 
paths for filament formation. The first - obvious - one 
is global collapse along the shorter axes. The second 
one - less obvious - arises from the fact that the radial 
accelerations in a two-dimensional elliptical or circular 
cloud of uniform density diverge at the edges, i.e. mate- 
rial at the edges experiences the strongest accelerations 
inwards, leading to a pil e-up o f gas at the edge, i.e to the 
formation of a filament (jBurkert fc Hartmannll2004f ). 

Obviously, "real" clouds are three-dimensional. How- 
ever, if clouds form in (laterally constrained) collid- 
ing flows, they will have a finite extent, and more im- 
portantly, global gravity will not have had sufficient 
time to lead to a centrally peaked density profile. The 
strong thermal instabilities lead to an initially rather 
thin sheet, only broadened by the dynamical instabili- 
ties. Thus, to zeroth order, in such a scenario the sheet- 
appr oximation is quite reasonable. T he numerical mod- 
els of'Vaz quez-Semadeni et all (|2007^ and the global col- 
lapse pattern of models Gfl and Gf2 above support this 
interpretation. In fact, the dynamical instabilities will 
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Fig. 10. — Top: One-dimensional velocity dispersions against 
time, along the inflow direction (ctx, thin lines), perpendicular to 
the inflows (tTyzmedium lines), and total (cr^{T < lOOK), thick 
lines), for all models. Bottom: One-dimensional velocity disper- 
sions against time. We distinguish between the velocity dispersion 
within each core, averaged over all cores {{dcore), thin lines), the 
velocity dispersion of the gas taken over all cores (dcora, medium 
lines), and the total velocity dispersion (as in top panel, thick lines). 



enhance a "crumpling" of the cloud in the lateral direc- 
tions once global gravity dominates. As the discussion in 
t )3 3.1 3.lTT] shows however, stars will have formed locally 
by then. 

It is reasonable to assume that the (idealized) inflows 



implemented in our models will not have a circular cross 
section. In galaxy mergers or in the collisions of super- 
shells in the LMC, the thickness of the disk would limit 
the thickness (vertical extent) of the flow. A similar as- 
sumption seems valid for spiral shocks in the Galaxy, 
where the disk pot ential would a gain lead to a flattening 
of the inflows (e.g. iDobbs et all [2006). 

4.2. Formation of Massive Cores 

Ripples in the flow collision interface can focus the in- 
streaming gas, leading to a very efficient mechanism to 
form massive cores (Figs. HI [S|). 10 Myr after flow col- 
lision, the flrst core in model Gs has a mass of approx- 
imately 150 M0 and a diameter of 1 pc, corresponding 
to a column density of 3.6 x 10^^ cm~^. Likewise, af- 
ter 10 Myr, the mean column density in the box along 
the inflow direction will be Ntot = 3 cm"'^ x (44 pc + 
7.9 km s"i X 10 Myr) = 1.1 x 10^^ cm"^. The dynam- 
ical focusing - helped along by cooling and eventually 
gravity - leads to an excess of column-density by a fac- 
tor of « 30. 

This is as good a place as any to remind the reader 
that we are leaving out the crucial step of molecule for- 
mation in our models. The flrst appearance of molecules 
sets the clock for the lifetime of the cloud. We can 
only mimic this by arguing that once column densities of 
« 10^ ^ cm~^ are reached, we regard the cloud as "molec- 
ular" . iHartmann et al.l (|2001h point out that this column 
density is of the same order as the critical column density 
for gravity to become dominant. 

With the NTSI as driving mechanism, the cores will 
form at a certain distance from the bulk of the cloud 
(Fig. 2]). This distance is basically given by the am- 
plitude the NTSI has reached up to that point. The 
somewhat peculiar location could have profound reper- 
cussions on the effect of feedback from the massive stars 
forming in the core: the winds and the expanding HII re- 
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Fig. 11. — Histograms of the Jeans length (left column) for all th ree self-g ravitating models, and scatter plots of the Jeans length against 
density, for all models at three times. The dashed lines denote the lTruelove ct al.. (.1997!) criterion with a safety factor of 2 and 4. 



gions or supcrnovac might lead to a further compression 
of the already massive cloud next to the young stars, thus 
triggering further star formation. In contrast, if the stars 
were located inside the bulk of the cloud, the stellar feed- 
back could be expected to provide partial support to the 
cl oud, or more likel y disperse the lower-density regions. 

iKeto et all ()2005( ) analyzed the molecular cloud popu- 
lation of the starburst galaxy M82, arguing that in many 
cases, star formation seems to proceed from the outside 
inwards, i.e. that the regions of massive star formation 
(indicated by HII regions and supernova remnants) are 
to be found at the edges of the molecular clouds. They 
argue that a sudden increase in external pressure (e.g. 
a shock wave traveling through the cloud) triggers star 
formation in a previously existing cloud. The dynam- 
ical focusing effect discussed here might offer another 
explanation, which then would imply that the central 
molecular cloud is still forming - an alternative which 
may be even more attractive in view of the notorious 



problem of stabilizing an object of many Jeans masses 
against gravitational collapse (iBurkert fc HartmannI 
l2ci04HBallesteros-Paredes fc Hartmannll2007D . 

The role of the "external pressure increase" in our mod- 
els is taken over by the dynamical focusing, i.e. by the ex- 
cess ram pressure at the troughs of the perturbed sheet. 
As Figure [T^] demonstrates, the thermal pressure does 
not vary strongly. 

Thus, the gas still evolves sub-isothermally i.e. the 
effective adiabatic exponent 7 < 1. 

At late times (top panel, 15.2 Myr), kinetic and grav- 
itational pressure are in equipartition. Once again, we 
emphasize that this does not mean that the system is in 
equilibrium: t he kinetic energy just foll ows the gravita- 
tional energy (|Ballesteros-Paredes|[2006l ) . 

4.3. Early Fragmentation and Core Mass Distributions 

The core mass distributions (Fig- P) are consistent 
with observations (jKramer et alj HOOGTlSchneider et al.l 
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20021). although th ey are flatter than the Salpeter IMF. 



Alves et all ()2007t ) quote a core mass distribution for the 
Pipe nebula close to the Salpeter IMF, however, as they 
point out, the cores that they are using for analysis are 
probably the direct progenitors of stars - a stage we can- 
not hope to reach with the models presented here. Bear- 
ing in mind that the distributions for models Gfl, Gf2 
and Gs arc probably slightly too flat for numerical rea- 
sons ( H'S 3.2 3.2.3p . their similarity to the non-gravitating 
model Hfl substantiates the claim that the density sub- 
structure - specifically the core mass distribution - in 
molecular clouds cou ld very well arise very early on dur- 
ing t heir formation (| Alves et al.l 120071 : iHennebelle et al.] 
I2007D - driven more by thermal than by gravitational 
fragmentation. 

However, we caution that we have not followed the 
cloud evolution through the stage of molecular hydrogen 



formation, an investigation we defer to a later paper. 
Finally, we again emphasize the difficulty of determin- 
ing mass spectra from simulations with limited dynamic 
range. 

4.4. Star Formation Duration and Stellar Age Spread 

From Figure \5\ one would infer an "age spread" of the 
stellar population in the complex of up to 8 Myr. This 
of course assumes that massive, self-gravitating cores will 
sit around for several Myr without forming massive stars 
whose energy input will tend to disrupt the cloud - a 
highly improbable supposition. While this aspect of our 
simulation is quantitatively unrealistic, it does suggest 
the following qualitative points: 

(1) Even in a completely dynamic simulation with- 
out magnetic or driven turbulent support, the overall 
timescales of star formation can be longer than the local 
collapse times. Putting this another way, the age spread 
in a star-forming region is an upper limit to the timescale 
of local collapse; it is only equal to the local collapse time 
if all star formation is globally synchronized, which is less 
and less plausible on larger and larger scales. 

(2) The large-scale flow picture of star-forming cloud 
accumulation is attractive in that it allows for the for- 
mation of stars over timescales (1-2 Myr) short com- 
pared with lateral cross ing times (10-20 Myr), as ob- 
served (jHartmann et al. 2001). This does not mean that 
a region of space may not have a signiflcant age spread. 
In the paradigm we are pursuing, where non-linear fluc- 
tuations are important, it is plausible that a few large 
perturbations might collapse long before more general 
star formation ensues, as we see in our simulations. Al- 
though we form mostly massive cores, other initial con- 
ditions - and higher numerical resolution (see 3.4p - 
may lead to the formation of a few low-mass stars ini- 
tially from the low-probability high density tail of pertur- 
bations, with the smaller, more frequent perturbations 
constituting the bulk of the star formation at a later 
time. Again, such age spreads provide no constraint on 
local collapse timescales. 

4.5. The Role of Magnetic Fields 

Magnetic fields are not included in the models pre- 
sent ed here. Based on tw o -dimen sional simulations 
by IVazguez-Semadeni et al] ()1995D . [Hartmann et al] 
(|200l[ ) envisaged the large-scale converging flows pref- 
erentially running alo ng (dynamically dominant, e.g. 
iHeiles fc Trolandl |2005[ ) magnetic field lines in the dif- 
fuse ISM, arguing that if the field lines were oriented 
perpendicularly, the stiffening of the equation of state 
due to the increased magnetic pressure would counter 
the lowering of the effective adiabatic index d ue to 
(atomic line) cooling (see also iBergin et al]|2004[ ). Al- 
though there are problems with this scenario (for one, 
substructures can form even in the case of perpendic- 
ular field lines and lead to further fragmentation [e.g. 
iHeitsch et al.l l2007| , and then, this scenario is clearly 
motivated by two-dimensional simulations, not allowing 
for interchange modes in the instabilities. These usu- 
ally grow at least at th e hydrodynamical rate, see e.g. 
IStone &: Gardiner! l2007f l , we will adopt it for the cur- 
rent study. Its consequence is that since material is be- 
ing piled up along the field lines, the mass-to-flux ratio 
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would e asily exceed the critical value f or magnetical dom- 
inance (|Mouschovias fc Spitzed [l976l ) after a few Myr, 
thus rendering the magnetic fields unimportant for the 
local gravitational collapse. It is true that the global 
dynamics of the cloud especially during its early evolu- 
tionary stages could be fundamentally influenced by the 
presence of magnetic fields. We will defer this important 
problem to a future paper. 

5. SUMMARY 

In an ext ension of our previous work (jHeitsch et al.l 
l2005l l2006f ) we presented three-dimensional models of 
molecular cloud formation in colliding flows including 
self-gravity. At an effective resolution of 512'^ they cur- 
rently are the most highly-resolved models of this kind. 
We used a fixed-grid code and a non-periodic Poisson- 
solver. The initial conditions emphasized the effects of 
global versus local gravity on the newly forming cloud. 

The model resolution did not allow us to follow the 
collapsed cores and to study their properties. Stellar 
feedback was not included. Thus, as time increases, 
the global evolution of the molecular cloud gets less and 
less realistic. Furthermore, we cannot make statements 
about the lifetime of the molecular cloud. Bearing these 
caveats in mind, we found that: 

(1) Any perturbation in the colliding flows is strongly 
amplified by a combination of thermal and dynamical 
instabilities. Specifically, the thermal effects lead to lo- 
cal high-density fragments that subsequently collapse be- 
fore the (still forming) cloud can collapse globally. Thus, 
cloud formation in colliding fiows allows the rapid on- 
set of local star formation while evading the problem of 
globally colla psing clouds. T his is consistent with earlier 
findings (Burk ert fc Hartma nn 2004) that local collapse 
can only win over global collapse in the presence of early 
non-linear density perturbations (Fig.[7])- 

(2) Even in the highly dynamical environment of the 
colliding flows, an elongated finite cloud under global 
gravity can form one or more filaments by collapsing 
along its shorter axes. This lateral collapse opens up 
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a further mass reservoir for star formation (Figs [l] [2|) . 

(3) Dynamical focusing, i.e. the deflection of incoming 
gas due to ripples in the interface between colliding gas 
streams, leads to very efficient high-mass core formation 
(Figs m El): core masses of a few 100 Mq are reached 
within « 10 Myr. 

(4) The clouds are not in any state of equilibrium at 
any time. Different scales submit to gravitational col- 
lapse at different times, with the small scales going first 
(Figs [SI [7]). Still, the global equipartition parameter 
(eq. [2]) ranges around 1 within a factor of a few - con- 
sistent with observations when taking into account the 
observational uncertainties. 

(5) The similarity of the core mass distribution for 
models with and without self-gravity and their robust- 
ness with time indicate that it might be set very early 
on during cloud formation. Thus, the core mass dis- 
tribution would be a consequence of fragmentation due 
to (magneto-)hydrodynamical instabilities and cooling 
rather than due to self-gravity (Fig. [9]). 

(6) The clouds become globally unstable in the absence 
of feedback, leading to an exploding "star formation ef- 
ficiency" (in quotes, since we only can form cores in our 
models, but not stars). The turbulence imparted by the 
(continuing inflows) does not suffice to balance global 
gravity. As mentioned in [JT] and above, the absence of 
feedback in our models does not allow us to test how 
long the clouds will survive. However, our models show 
that the clouds will not disperse "on their own accord" 
in the background flows, i.e without the help of stellar 
feedback. 
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